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Abstract 

An ubiquitous property of biological sensory systems is adaptation: a step increase in 
stimulus triggers an initial change in a biochemical or physiological response, followed by 
a more gradual relaxation toward a basal, pre-stimulus level. Adaptation helps maintain 
essential variables within acceptable bounds and allows organisms to readjust themselves 
to an optimum and non-saturating sensitivity range when faced with a prolonged change 
in their environment. Recently it was shown theoretically and experimentally that many 
adapting systems, both at the organism and single-cell level, enjoy a remarkable additional 
feature: scale invariance, meaning that the initial, transient behavior remains (approxi- 
mately) the same even when the background signal level is scaled. In this work, we set 
out to investigate under what conditions a broadly used model of biochemical enzymatic 
networks will exhibit scale- invariant behavior. An exhaustive computational study led us 
to discover a new property of surprising simplicity and generality uniform linearizations 
with fast output (ULFO), whose validity we show is both necessary and sufficient for scale 
invariance of enzymatic networks. Based on this study we go on to develop a mathematical 
explanation of how ULFO results in scale invariance. Our work provides a surprisingly 
consistent, simple, and general framework for understanding this phenomenon, and results 
in concrete experimental predictions. 

Author summary: An ubiquitous property of biological sensory systems is adaptation: a step 
increase in stimulus triggers an initial change in a biochemical or physiological response, followed 
by a more gradual relaxation toward a basal, pre-stimulus level. Adaptation helps maintain 
essential variables within acceptable bounds and allows organisms to readjust themselves to 
an optimum and non-saturating sensitivity range when faced with a prolonged change in their 
environment. Recently, it was shown theoretically and experimentally that many adapting 
systems, both at the organism and single-cell level, enjoy a remarkable additional feature: 
scale invariance, meaning that the initial, transient behavior remains (approximately) the same 
even when the background signal level is scaled. In this work, we develop a mathematical 
characterization of biochemical enzymatic networks that exhibit scale-invariant behavior and 
make concrete experimental predictions. 
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1 Introduction 



The survival of organisms depends critically upon their capacity to formulate appropriate re- 
sponses to sensed chemical and physical environmental cues. These responses manifest them- 
selves at multiple levels, from human sight, hearing, taste, touch, and smell, to individual 
cells in which signal transduction and gene regulatory networks mediate the processing of mea- 
sured external chemical concentrations and physical conditions, such as ligand concentrations 
or stresses, eventually leading to regulatory changes in metabolism and gene expression. 

An ubiquitous property of biological sensory systems at all levels is that of adaptation: a step 
increase in stimulus triggers an initial, and often rapid, change in a biochemical or physiological 
response, followed by a more gradual relaxation toward a basal, pre-stimulus level [I]. Adapta- 
tion plays a role in ensuring that essential variables stay within acceptable bounds, and it also 
allows organisms to readjust themselves to an optimum and non-saturating sensitivity range 
even when faced with a prolonged change in their operating environment, thus making them 
capable of detecting changes in signals while ignoring background information. 

Physiological examples of adaptation in higher organisms include phenomena such as the control 
of the amount of light entering eyes through the contraction and relaxation of the pupil by the 
nervous system, which brings intensities of illumination within the retinal working range, or the 
regulation of key metabolites in the face of environmental variations [20]. At the single-cell level, 
one of the best understood examples of adaptation is exhibited by the E. coli chemotaxis sensory 
system, which responds to gradients of nutrient and ignores constant (and thus uninformative) 
concentrations [39] . The term "exact" or "perfect" adaptation is employed to describe 
processes which, after a transient, return with very high accuracy to the same input-independent 
level. In practice, however, an approximate adaptation property is usually adequate for proper 
physiological response [30J. 

By definition, neither the concepts of perfect nor approximate adaptation address the charac- 
teristics of the transient signaling which occurs prior to a return to steady state. The amplitude 
and other characteristics of transient behaviors, however, are physiologically relevant. In this 
more general context, a remarkable phenomenon exhibited by several human and animal sen- 
sory systems is scale invariance or logarithmic sensing [20] [22] [48]. This means that responses 
are functions of upon ratios (in contrast to actual magnitudes), of a stimulus relative to the 
background. There is evidence for this phenomenon at an intracellular level as well. It appears 
in bacterial chemotaxis [18] [35], in the sensitivity of S. cerevisiae to fractional rather than 
absolute pheromone gradients [31], and in two mammalian signaling systems: transcriptional 
as well as embryonic phenotype responses to /3-catenin levels in Wnt signaling pathways [T3] , 
and nuclear ERK localization in response to EGF signaling [10]. Scale invariance allows sys- 
tems to react to inputs ranging over several orders of magnitude, and is speculated to help 
make behaviors robust to external noise as well as to stochastic variations in total expressed 
concentrations of signaling proteins [41] . 

Mathematically, scale invariance is defined by the following property of transient behaviors 
[H]: if a stimulus changes from a background level uq to a new level u, then the entire time 
response of the system is the same as if the stimulus had changed, instead, from a background 
level puo to pu. In other words, only the ratio (or "fold-change") pu/puo = u/uq is relevant 
to the response; the "scale" p is irrelevant. For this reason, the term "fold change detection" 
is interchangeably used instead of scale- invariance. Scale invariance implies adaptation, but 
not every adaptive system is scale invariant [H]. A mathematical analysis of scale-invariance 
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was initiated in [JT], |40j . Predictions regarding scale-invariance of E. coli chemotaxis were 
subsequently experimentally verified [23]. While adaptation can be often understood in terms 
of control-theoretic tools based on linearizations [32] [52] [35] [E] [25], scale invariance is a 
genuinely nonlinear property; as a matter of fact, a linear system can never display scale- 
invariance, since the response to an input scaled by p will also be scaled by this same factor 
p. 

In this work, we focus on enzymatic signal transduction systems, which involve the activa- 
tion/deactivation cycles that typically mediate transmission of external signals to transcription 
factors and other effectors. Networks involving such enzymatic cycles are involved in signal 
transduction networks from bacterial two-component systems and phosphorelays jSJIH] to actin 
treadmilling [9], guanosine triphosphatase cycles [IT], glucose mobilization [19], metabolic con- 
trol [15], cell division and apoptosis [l6j, cell-cycle checkpoint control [24J, and the eukaryotic 
Mitogen- Activated Protein Kinase (MAPK) cascades which mediate growth factor inputs and 
determine proliferation, differentiation, and apoptosis [^ IHf fT5l 150} 15]. 

Given the biological importance of these processes, and the already observed scale-invariance in 
some of these pathways [T5] [10] , we pose here the following question: which enzymatic networks 
do not merely adapt, but also display scale invariance? In order to answer this question, 
we performed an exhaustive computational study of all 3-node networks, finely sampled in 
parameter space. Only about 0.01% of these networks are capable of (approximate) adaptation. 
Testing which of these adapting networks also display scale- invariant behavior, we found that 
only about 0.15% of them did. Once that this small subclass was identified, we turned to 
the problem of determining what network characteristics would explain the results of these 
numerical experiments. We discovered a surprisingly simple and general property, which we 
call uniform linearizations with fast output (ULFO), that is displayed by all the networks in 
this subclass, and here we provide a theoretical framework that explains conceptually why this 
property is both necessary and sufficient for scale invariance of enzymatic networks. As an 
application, we consider a recently published model [17] of an eukaryotic enzymatic system, 
specifically the pathway involved in the social amoeba Dictyostelium discoideum's chemotactic 
response to cAMP. and show that our conditions are satisfied in appropriate ranges of cAMP 
input. 

Characterizations of this sort allow one to understand which networks are robust to scale 
uncertainty, and constitute a powerful tool in allowing one to discard putative mechanisms that 
are not consistent with experimentally observed scale-invariant behaviors [40j, [23] . 

2 Results 

2.1 Three-node enzymatic networks 

We consider networks consisting of three types of enzymes, denoted respectively as A, £>, and 
C. Each of these enzymes can be in one of two states, active or inactive. The fractional 
concentration of active enzyme A is represented by a variable xa — so xa — 1 — xa is 

the fraction of inactive enzyme A. Similar notations are used for B and C. Only enzyme A is 
directly activated by an external input signal, and the response of the network is reported by the 
fraction of active C. Enzyme B acts as an auxiliary element. Each enzyme may potentially act 
upon each other through activation (positive regulation), deactivation (negative regulation), or 
not at all. If a given enzyme is not deactivated by any of the remaining two, we assume that it is 
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constitutively deactivated by a specific enzyme; similarly, if a given enzyme is not activated by 
any other, there is a constitutively activating enzyme for it. One represents networks by 3-node 
directed graphs, with nodes labeled A, £>, C, and with edges between two nodes labeled + and 
— (or "— and "H") to denote positive or negative regulation respectively; no edge is drawn if 
there is no action. There are 3 2 = 9 potential directed edges among the three nodes (A to A, A 
to £>, etc.), each of whose labels may be +, — , or "none" if there is no edge. This gives a total 
of 3 9 = 19, 683 possible graphs. One calls each of these possible graphs a topology. Discarding 
the 3,645 topologies that have no direct or indirect links from the input to the output, there 
remain 16,038 topologies. 

2.2 Specification of a dynamic model 

We quantify the effects of each existing regulatory interaction by a Michaelis-Menten term and 
write a three- variable ordinary differential equation (ODE) that describes the time evolution of 

XA(t), XB(t), and xc(t): 



The ICs denote Michaelis-Menten, and the fc's catalytic, rate constants associated to each 
regulatory interaction. All the summations range over i = 1,...,6. Each "T^" represents 
one of A, £>, C, Ea, Eb, Eq, the activating enzymes in the respective equations, and each 
"Wi" one of A, £>, C, Fa, Fb, Fq, the deactivating enzymes; E and F are the constitutively 
activating and deactivation enzymes, buffered at constant concentrations. (Lower-case variables 
Vi,W{ = xa, • • • ,%F C denote active fractions) As an exception, the equation for node A does 
not include an Ea term, but instead includes a term kjjAU - — that models activation of A 



by an external input whose strength at time t is given by u — u(t) and whose values u(t) stay 
within a range [m, No enzyme appears both an activator and as a deactivator of any given 
component, that is, kx^ky^ — 0, kx^ky^ — 0, and kx^ky^ — 0, and constitutive enzymes 
are included only if the reaction would be otherwise irreversible. For example, the topology 
shown in Fig{T] is described by the following following set of ODE's: 
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Figure 1: Topology 2293 
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k UA U • X A k B AXB • XA kcAXC 'XA , n , 

xa — ~ (2aJ 

x A + K UA x A + K B A x A + K C a 

k A BXA • x B k FB BX FB • x B 

xb — — 77 — V^bj 

xb + Kab x b + Kf b b 

kACXA • Xc k B C • X B XQ kccXC ' %C , x 

X(j — — — (Ac) 

xc + K A c xc + K B c xc + K C c 

The term circuit is used to refer to a given topology together with a particular choice of the K 
and k parameters. The three-node model in Eq|l]was employed by Ma et al. [25J, in order to 
classify the minimal enzymatic circuits that adapt. (With the model in [25 j that we adopted, 
there is no direct connection from the input to the output node, and two-node networks are 
not sufficient for adaptation, while larger adapting networks contain these three-node networks 
[25] . If one allows direct connections from input to outputs, then two-node networks are able 
to display adaptation.) The same paradigm has since been used to investigate other network 
characteristics as well [38], [51] . 



2.3 Adaptation 

Following [12] , we define adaptation behavior in terms of two functional metrics. The first metric 
quantifies the following effect: if we start at steady state, and then step the input at time t — 
from a value uo to a different constant value u\, then the system's output, as reported by a 
response variable y{t) (where y(t) — xc(t) in Eq{T]), should return asymptotically to a value 
that is close to the original value y(0). The relative difference in initial and final response 
A^° = \y(+oo) — 2/(0) | provides a measure of adaptation precision. We say that a system 
is (approximately) adaptive provided that, for all inputs in the valid range, A^/A u < 0.1, 
where A u = \u± — uo\ / \uo\ is the relative change in input. In particular, exact or perfect 
adaptation means that = 0. The 10% error tolerance is natural in applications, and 
the qualitative conclusions are not changed by picking a smaller cutoff [25J. A second metric 
relies upon the maximal transient difference in output, normalized by the steady- state output, 
A™ ax = max |?/(£) — y(0)| / |y(0)|. A signal- detection property for adaptation [33], [2], should 
be imposed in order to rule out the trivial situation A™ ax ^ in which a system's output is 
independent of the input. To avoid having to pick an arbitrary threshold, in this study we 
follow the convention in [25] of requiring the sensitivity A™ ax /A u to be greater than one. 



2.4 Scale invar iance 



Scale invariance is the property that if a system starts from a steady state that was pre- 
adapted (t < 0) to a certain background level i^o, and the input is subsequently set to a 
new level u at t = 0, then the entire time response of the system y UQ ,u{t) is the same as the 
response y P u ,pu(t) that would result if the stimulus had changed, instead, from puo to pu. This 
property should hold for scale changes p > that respect the bounds u < u < u on inputs. 
For example, recent microfluidics and FRET experimental work [23J verified scale-invariance 
predictions that had been made in [41] for bacterial chemotaxis under the nonmetabolizable 
attractant o:-methylaspartate (MeAsp) as an input. In these experiments, E. coli bacteria 
were pre-adapted to input concentrations and then tested in new nutrient gradients, and it was 
found experimentally that there were two different ranges of inputs [?/ l5 Hi] and [u 2l U2] m which 
scale-invariance holds, the "FCD1" and "FCD2" regimes, repectively. (The term fold-change 
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detection, or FCD, is used to reflect the fact that only the ratio or fold-change pu/puo = u/uo 
can be detected by the response y(t).) More generally, the mathematical definition of (perfect) 
scale invariance [40J imposes the ideal requirement that the same response invariance property 
is exhibited if u — u(t), t > is any time- varying input. The experiments in [23] included 
excitation by certain oscillatory inputs, for example. In practice, however, this property will 
always break down for high-frequency inputs, since there are limits to the speed of response of 
biological systems. 

2.5 Adaptive systems need not be scale-invariant 

As an illustration of a (perfectly) adaptive yet not scale-invariant system, consider the following 
equations: 

xa — k\u — k2XB (3a) 
xb — ksXA — k^XB (3b) 
x c = hx A - k 6 x B x c (3c) 

which is a limiting case of the system described by Eqj2]when kcA, kcci Kjja, Kb a, Kab, Kac ~ 
0, k B c = k 6 K B c, K B c > 1 (so -k B cXBXc/(xc+K B c) ~ -k 6 x B x c ), and k FBB x FB = k 2 K FBB 
and K Fb b ^> 1. This network perfectly adapts, since at steady state the output is xc — xc — 
^4^5/(^3^6)5 no matter what is the magnitude of the constant input u, and in fact the system 
returns to steady state after a step change in input u, with xc(t) —> xc as t — >> 00 (general sta- 
bility properties of feedforward circuits shown in [33] )• On the other hand, the example in Eqj3] 
does not display scale invariance. Indeed, consider the solution from an initial state pre-adapted 
to an input level uo, that is xa(0) = ^1^4^0/(^2^3)5 xb(0) — fci^o/^ and xc(0) = ^4^5/(^3^), 
and the input u(t) = u\ for t > 0. Then, xc(t) = k^/^k^ko) + k\k^{ui — uo)t 2 /2 + 0(t 3 ) 
for small t > 0. Since the t 2 coefficient in this Taylor expansion gets multiplied by p when 
uo is replaced by puo and u\ is replaced by pu\, it follows that the transient behavior of the 
output xc(t) depends on p. Interestingly, if the equation for the third node is replaced by 
xc — k$XA/xB — kftXc, that is to say the activation of C is repressed by A, instead of its 
de- activation being enhanced by A, then scale invariance does hold true, because x^(t) and 
xs(t) both scale by p when uq 1— )► puo, u\ 1 — ^ 1x0, and C{t) depends on the ratio of these two 
functions (in particular, the t 2 /2 term is ^2^5(1^1 — uq)/uq). Such a repression is typical of 
genetic interaction networks, but is not natural in enzymatic reactions. 

It turns out that the example described by Eqj3]is typical: no enzymatic network described by 
Eq{T]can display perfect scale-invariant behavior. This fact is a consequence of the equivariance 
theorem proved in [ID] (see Materials and Methods). Thus, a meaningful study of enzymatic 
networks, even for perfectly adaptive ones, must rely upon a test of approximate scale invari- 
ance. Instead of asking that yuo,u(t) — Upuo,pu 

(£), as was the case in the theory developed in 
[H] [30], one should require only that the difference be small. To investigate this issue, we 
computationally screened all 3-node topologies through a high-throughput random parameter 
scan, testing for small differences in responses to scaled steps. We found that approximately 
0.01% of the samples showed adaptation, but of them, only about 0.15% passed the additional 
criterion of approximate scale invariance (see Materials and Methods). These samples belonged 
to 21 (out of 16,038 possible) topologies. As an example of the behavior of one of these, Fig(2] 
shows a response resulting from a 20% step, from 3 to 3.6, compared to the response obtained 
when stepping from 5 to 6; the graphs are almost indistinguishable. (See SI Text for an enu- 
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meration of circuits and corresponding plots) . In the following discussion, we will refer to these 
surviving circuits, and their topologies, as being "approximately scale invariant" (ASI). 

We found that all ASI networks possess a feedforward motif, meaning that there are connections 
A —> B —> C and as well as A —> C. Such feedforward motifs have been the subject of 
extensive analysis in the systems biology literature [lj. and are often involved in detecting 
changes in signals [28]. They appear in pathways as varied as E. coli carbohydrate uptake 
via the carbohydrate phosphotransferase system [2Tj . control mechanisms in mammalian cells 
[26] . nitric oxide to NF-^B activation [271 EH], EGF to ERK activation [33 [32], glucose to 
insulin release [311 EH], ATP to intracellular calcium release [36J, and microRNA regulation 
[49]. The feedforward motifs in all ASI networks are incoherent, meaning such that the direct 
effect A C has an opposite sign to the net indirect effect through B. An example of an 
incoherent feedforward connection is provided by the simple system described by EqjSJ where 
the direct effect of A on C is positive, but the indirect effect is negative: A activates B which 
in turn deactivates C. (Not every incoherent feedforward network provides scale invariance; a 
classification of those that provide exact scale invariance is known [40].) It is noteworthy that 
all ASI circuits have a positive regulation from A to B and a negative regulation from B to A. 

We then discovered a surprising common feature among all ASI circuits. This feature can best 
be explained by a further examination of the example in Eqj3j 
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Figure 2: Scale-invariance: plots overlap, for responses to steps 3—^1.2*3 and 5—^1.2*5. Network 
is the one described by Eqj2j Random parameter set: K UA =om39i8 k UA =iiA472i9, K BA =o.ooi688 

/c£> A =44.802268, K C A=90.209027 /c C a=96.671843, X A£ >=0.001191 fc Afl =1.466561, Kf b =9.424319 kp B =22.745736, 
K AC =0.113697 fc AC =l-211993, 1^=0.009891 /c B c=7.239357, X C c=0.189125 fc C c=l 7.910182 



2.6 Approximate scale invariance 

Continuing with example in Eqj3j let us suppose that ^1,^2,^3,^4 <C ^5,^65 so that the output 
variable y — xc reaches its steady state much faster than xa and xb do. Then, we may 
approximate the original system by the planar linear system represented by the differential 
equations for xa and xb together with the new output variable y(t) = h(xA(t),XB(t)) = 
kxA(t)/xB(t), where k = k^/k^. This reduced planar system, obtained by a quasi-steady state 
approximation, has a perfect scale-invariance property: replacing the input u by pu results in the 
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solution (pxA(t),pxB(t)), and thus the output is the same: h(xA(t),XB(t)) = h(pxA(t),pxB(t)). 
The exact invariance of the reduced system translates into an approximate scale invariance 
property for the original three-dimensional system because, except for a short boundary-layer 
behavior (the relatively short time for xq to reach equilibrium), the outputs of both systems 
are essentially the same, y(t) w y(t). 



2.7 Generality of the planar reduction 

We found that, just as in the example in Eq(3]when ^1,^2,^3,^4 <C &5,&6, m ever y ASI cir- 
cuits the time scale of node C is much shorter than that of A and B. Therefore, the same 
two-dimensional reduction is always valid. It follows that one can drop the last equation, ap- 
proximating these circuits by planar systems that are described by only the two state variables 
xa and xb, where every occurence of xc in the first two equations of the right-hand side of 
Eq|I]is replaced by 1i(xa,xb), the function obtained by setting the right-hand side of the third 
equation in Eq|I]to zero and solving for the unique root in the interval [0, 1] of the quadratic 
equation. This reduced system, with y(t) = h(xA(t),XB(i)) as an output, provides an excellent 
approximation of the original dynamics. Figj3] compares the true response with the response 
obtained by the quasi-steady state approximation, for one ASI circuit (see SI Text for all com- 
parisons) . 
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C 6 (t) by solving quadratic eqn 

C 6 (t) as a solution of nonlinear diff. eqn. 
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Figure 3: QSS quadratic approximation. Network is the one described by Eq(2j Random 
parameter set is as in Fig|2| 



2.8 Generality of dependence on x A /x B 

In the example given by Eqj3j there were two additional key mathematical properties that 
made the planar reduction scale-invariant (and hence the original system approximately so). 
The first property was that, at equilibrium, the variable xq must be a function of the ratio 
%a/%b, and the second one was that each of xa and xb must scale by the same factor when 
the input scales by p. Neither of these two properties need to hold, even approximately, for 
general networks. Surprisingly, however, we discovered that both are valid with very high 
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accuracy for every ASI circuit. The equilibrium value of xc is obtained from setting the last 
right-hand side of Eq|I]to zero and solving for xc- A solution xc — h(xA,xs) in the interval 
[0, 1] always exists, because at xc = one has xq — 1 and thus the term is positive, and at 
xc = 1 one has xc = and so the term is negative. This right-hand side has the general form 
xa<I>(xc) + xb^(xc) + k(xc, xe c , xp c ), where (j) and 7 are increasing functions, each a constant 
multiple of a function of the form xc/{xc + K) or —xc/(xc + K). If the term n is negligible, 
then xa^(xc) + xb^{xc) — means that also (xa/xb)(/>(xc) +7(^c) = 0, and therefore xc at 
equilibrium is a (generally nonlinear) function of the ratio xa/xb- There is no a priori reason 
for the term k to be negligible. However, we discovered that in every ASI circuit, n w 0. More 
precisely, there is no dependence on the constitutive enzymes, and this "self-loop" link, when 
it exists, contributes to the derivative xc much less than the xa and xb terms, see FigjiJ 
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Figure 4: Relative contribution of terms in the equation for node C. The first two terms range 
in [—0.25,0.25] but self-loop magnitude is always less than 10 -3 . i.e. contribution or self-loop 
to xc is less than 1%. Similar results hold for all ASI circuits. Network is the one described by 
Eqj2j Random parameter set is as in Figj2j Similar results are available for all ASI circuits. 



2.9 Generality of homogeneity of x A ,x B 

The last ingredient of the example given by Eqj3]that plays a role in approximate scale invariance 
is that each of xa and xb must scale proportionately when the input is scaled. In that example, 
the property holds simply because the equations for these two variables are linear. In general, 
however, the dynamics of (xa, xb) are described by nonlinear equations. Thus it is remarkable 
that, in all ASI circuits, the property holds. We tested the property by plotting x Ait) /x sit) in 
a set of experiments in which a system was pre-adapted to an input value uq and the input was 
subsequently set to a new level u at t = 0. When going from puo to pu, we found that the new 
value x Ait) /x sit) was almost the same, meaning that xa and xb scaled in the same fashion. 
A representative plot is shown in Figj5j 



2.10 A new property: uniform linearizations with fast output 

The (approximate) independence of x Ait) /x sit) on input scalings is not due to linearity of 
the differential equations for xa and xsit). Instead, the analysis of this question led us to 
postulate a new property, which we call uniform linearizations with fast output (ULFO). To 
define this property, we again drop the last equation, and approximate circuits by the planar 
system that has only the state variables xa and xb, where every occurence of xc in their 
differential equations shown in Eq|l] is replaced by h{xA>>XB)- We denote by f(xA,xs,u) = 
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Figure 5: Constant A/B ratio in responses to 3—^1.2 * 3 and 5—>1.2 * 5. Network is the one 
described by Eqj2j Random parameter set is as in Figj2j Similar results are available for all 
ASI circuits (see SI Text). 



(fi(xA, xb, u), f2(%A, u)) the result of these substitutions, so that the reduced system is 
described in vector form by x = f(x,u), x — (xa,xb)- We denote by a(u) the unique 
steady state corresponding to a constant input u, that is, the solution of the algebraic equation 
f(a(u),u) = 0. We denote by A(u) = (df /dx)(a(u)) the Jacobian matrix of / with respect 
to x, and by B{u) — (df/du)(a(u)) the Jacobian vector of / with respect to u. The property 
ULFO is then defined by requiring time-scale separation for xc, that Ji(xa,xb) depends only 
on the ratio xa/xb, and: 



a(pu) = p<r(u), A(u) = A(v), B(u) = B(v) 



(4) 



for every u , v, and p such that u, v, and pu are in the range [u,u\. Notice that we are not 
imposing the far stronger property that the Jacobian matrices should be constant. We are only 
requiring the same matrix at every steady state. The first condition in Eqj4] means that the 
vector (j{u)/u should be constant. We verified that this requirement holds with very high accu- 
racy in every one of the ASI circuits. With u — 0.3 and u — 0.6, we have the following a(u)/u 
values, rounded to 3 decimal digits: (0.195,0.239), (0.193,0.237), (0.192,0.236), (0.191,0.235) 
when u — 0.3, 0.4, 0.5, and 0.6 respectively, for the network described by Eqj2]and the random 
parameter set in Figj2j Similar results are available for all ASI circuits (see SI Text). The Jaco- 
bian requirements are also verified with high accuracy for all the ASI circuits. We illustrate this 
with the same network and parameter set. Let us we compute the linearizations Aqs — ^4(0.3), 
Aqa — A(0A), . . . , £>o.6 = 23(0.6) and the average relative differences 



err 

ij 



E 

w=0.3,0.4,0.5,0.i 



(A u )ij - (^4o.45)ij 



(AoAb)ij 



and we define similarly £> err . These relative differences are very small (shown to 3 decimal 
digits) : 

ierr _ /^0.069 0.004\ err _ /0.002 N 
0.005 ) ' I 



A e 



thus justifying the claim that the Jacobians are practically constant. Similar results are available 
for all ASI circuits (see SI Text). 
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The key theoretical fact is that the property ULFO implies approximate scale-invariance, see 
Materials and Methods. 

2.11 A concrete example 

In a recent paper [17] Takeda and collaborators studied the adaptation kinetics of a eukary- 
otic chemotaxis signaling pathway, employing a microfluidic device to expose Dictyostelium 
discoideum to changes in chemoeffector cyclic adenosine monophosphate (cAMP). Specifically, 
they focused on the dynamics of activated Ras (Ras-GTP), which was in turn reported by RBD- 
GFP (the Ras binding domain of fluorescently tagged human Rafl), and showed almost perfect 
adaptation of previously unstimulated cells to cAMP concentrations ranging from 10 -2 nM to 
1 fiM. Furthermore, inspired by [25], the authors proposed alternative models for adaptation, 
and concluded that the best fit was obtained by using an incoherent feedforward structure. The 
model that they identified is given by the following system of 6 differential equations: 

k Rl (v + ri)(^ ot - R{) - k- Rl R! 
k R2 (v + r 2 )(4 ot - R 2 ) - k_ R2 R 2 

R\ + R 2 

k G EF u - k-GEFGEF 

k GA p u - kcApGAP 

k RAS GEF (RAS tot - Ras GTP ) - k- RAS GAP Ras GTP 
k°^ D (RBD tot - RBD^) - k™ BD Ras GTP RBD^ . 

The symbol v stands for the chemoeffector cAMP, and the authors assumed the existence of two 
different receptor populations (R\ and i?2, with very different K^s) which when bound pool 
their signals to downstream components (through u). The constants r\ and r 2 represent levels 
of constitutive activation. The variables GEF and GAP represent activation and deactivation 
of RasGEF and RasGAP, Ras GTP represents the activated Ras, and RBD C ^ describes the 
cytosolic reporter molecule RBD-GFP. Fig. [6] shows a schematic of the main players. 




Figure 6: The system studied in [47] 

The best-fit parameters obtained in [47J are as follows: R\ ot =0.1, i?2 0t = 0.9, r\ = 0.012nM, 
r 2 = 0.115nM, k Rl = 0.00267nM" 1 sec~ 1 , k- Rl = 0.16sec-\ k R2 = 0.00244nM" 1 sec~ 1 , k- R2 = 
l.lsec -1 , k GE F = 0.04sec _1 , k- G EF = 0.4sec _1 , k GA p = O.Olsec -1 , k GA p = O.lsec -1 , RAS tot = 



dRi 

~df 
dR 2 

~df 

u 

dGEF 

dt 
dGAP 

dt 

dRas GTP 
dt 

dRBD^ 
Jt 
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1, kjus = 390sec-\ k-RAS = 3126SCC- 1 , RBD tot = 1, k™ D = 0.53sec-\ k™ D = l.Osec" 1 . 
With these parameters, and cAMP concentrations which are small yet also satisfy r\ <C v(t) 
and T2 <C v(t), it follows that R\ « kR 1 R\°^v — k-R 1 R± and R2 ~ kR 2 R^^v — k-R 2 R2, so we 
may view u(t) as an input (linearly dependent on the external v(t)) to the three- variable system 
described by xa = GEF, xr = GAP, xc = Ras GTP . Since RBD C ^ depends only on xc, we 
may view xc as the output. This three- variable system (interpreted as having limiting values of 
Michaelis-Menten constants) has the ULFO property provided that the dynamics of xc are fast 
compared to xa and xr, which the identified parameters insure. So, we expect scale-invariant 
behavior. Indeed, Figj7| shows a simulation of the entire six-dimensional system (not merely 
of our 3-dimensional reduction) when using a step from 1 to 2 nM of cAMP, and shows that 
essentially the same response is obtained when stepping from 2 to 4 nM. This prediction of 

0.635 
U 8.5 
0.625 
82 
I 0.815 
0.81 
0.605 
0.8 
0.795 

10 20 30 40 50 60 70 60 90 100 

Figure 7: Scale-invariance for model from [47J: responses to steps 1— >2 and 2—^4 coincide 
scale-invariant behavior is yet to be tested experimentally. 

3 Discussion 

Work in molecular systems biology seeks to unravel the basic dynamic processes, feedback 
control loops, and signal processing mechanisms in single cells and entire organisms, both for 
basic scientific understanding and for guiding drug design. One of the key questions is: how 
can one relate phenotype (function) to interaction maps (gene networks, protein graphs, and 
so forth) derived from experimentation, especially those obtained from high-throughput tools? 
Answers to this question provide powerful tools for guiding the reverse-engineering of networks, 
by focusing on mechanisms that are consistent with experimentally observed behaviors, and, 
conversely, from a synthesis viewpoint, allow one to design artificial biological systems that are 
capable of adaptation [6] and other objectives. In particular, scale-invariance, a property that 
has been observed in various systems [13], [TO], can play a key role in this context, helping to dis- 
card putative mechanisms that are not consistent with experimentally observed scale-invariant 
behaviors [23]. Through a computational study, we identified a set of simple mathematical 
conditions that are used to characterize scale invariant enzymatic networks. 



■■■Input scaling p=1 
Input scaling p=2 
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4 Materials and Methods 



4.1 Computational screen 

We generalized and extended the computational protocol developed for adaptation in [25] to 
an investigation of approximate scale invariance. MATLAB® scripts were used, in conjunction 
with the software developed in [25]. In order to test inputs in ranges of the form a < u(t) < 2a, 
redefining the constant kjjA if needed, we take simply u — 0.3 and u — 0.6. We considered 
160,380,000 circuits, obtained from the 16,038 nontrivial 3-node topologies, each one with 10,000 
parameters sampled in logarithmic scale using the Latin hypercube method [17J. (We picked 
the ranges /c ca t=0.1-10 and K m =0. 001-100. A finer sampling does not affect conclusions in any 
significant way [25].) Of these, 0.01% (16,304) circuits showed adaptation, meaning that, as in 
[25] , when making a 20% step from uq = 0.5 to u\ — 0.6 the precision is 10% or better, and the 
sensitivity is at least unity. Approximate scale invariance (ASI) was then tested by also perform- 
ing a 20% step experiment from uq = 0.3 to u\ — 0.36 and requiring that the relative difference 
between the responses be at most 10%: max^ {| 2/0.6 (£) — 2/3.6 (^) I / m & x (2/o.6(£) — 2/3.6 (£))} < 0.1 
Of the adapting circuits, about 0.15% (25 circuits, classified into 21 different topologies) were 
determined to be ASI. 

4.2 ULFO implies approximate scale invariance 

Consider a system of n differential equations with input signal u, 

x — f(x, u) 

with the variables x evolving on some closed bounded set and / differentiate, and suppose 
that for each constant input u* there is a unique steady state x* = a(u*) with the conditions 
in Eqj4] and an output 

y(t) = h(x(t)) 

such that h is differentiate and homogeneous of degree zero (h(px) = h(x) for nonzero p). We 
view 3-node enzymatic networks as obtained from a set of n + 1 equations 

x — F{x,z,u) 
ez — G(x,z) 

with n — 2, x — (xa,xb), an d z = xq (0 < e <C 1 represents the faster time scale for xc), and 
we are studying the reduced system x = f(x, u) — F(x, a(x),u) obtained by solving G(x, z) = 
for z — ol{x) and substituting in F. Consider a time interval [0,T], a constant input u*, and 
a possibly time- varying input u(t), t > 0, as well as a scaling p > 0, such that all values u*, 
pu*, u(t), pu(t) are in the input range [u,u]. The solutions of x = f(x,u) with initial condition 
x(0) = cr(u*) and of i = f(z,pu) with initial condition z(0) = a(pu*) are denoted respectively 
by x(t) and z(t), and the respective outputs are y(t) = h(x(t)) and y p (t) = h(z(t)). We wish to 
show that these two responses are approximately equal on < t < T. Write 5(t) = u{t) — u*. 
From Theorem 1 in [12] we know that 

x(t)=x(0) + Z(t) + o(\\6\\) 
where \\S\\ = sup 0<t<T \S(t)\ and £ is the solution of the variational system 

i(t) = A£(t)+B5(t) 
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with £(0) = 0, and that 

z(t) = z(0) + C(t) + o(\\pu - pu*\\) = z(0) + C(t) + o(\\S\\), 

where 

C = Aat) + Bp8{t) 

with C(0) = 0. By linearity, £ = Using z(0) = a(pu*) = pa(u*) = px(0), we have that 
px(t) - z(t) = o(\\6\\). Thus, 

y(t) = h(x(t)) = h(px(t)) = h(z(t) + o(\\6\\)) . 

If K is an upper bound on the gradient of /i, then 

\y p (t) - y(t)\ = \h(z(t)) - h(z(t) + o(||*||))| < Ko(\\6\\)- 

Thus, the relative error sup t \y p (t) — y(t) \ / sup t \u(t) — u* \ converges to zero as a function of the 
input perturbation u(t) — u*. As a numerical illustration, we consider again the the network 
described by Eqj2] and the random parameter set in Figj2j We compare the relative error 
between the original nonlinear system, with initial state £ = (xa,xb) corresponding to u — 0.3, 
and applied input u — 0.36, and the approximation is £ + where the z solves the linear 
system with initial condition zero and constant input 0.06. The maximum approximation error 
is about 5% (to 3 decimal places, 0.055 for xa and 0.01 for xb)- When stepping from u = 0.5 
to u — 0.6, the error is less than 3% (0.028 and 0.005 respectively). Similar results are available 
for all ASI circuits (see SI Text). 

4.3 Impossibility of perfect scale- invar iance 

Consider any system with state x — (xa,xb, xq), output xc, and equations of the general form 
xa = f(x) + G(xa)u, x b = g(x), x c = h(x) = x A a(x c ) + x B b(x c ) + c(x c ). 

x A = f(x) + G(x A )u 
x B g(x) 

xc = h(x) = xacl(xc) + x B b(xc) + c(xc) • 

It is assumed that a(xc) ^ for all xc, G(xa) ^ for all xa, G := sup x G(x) < oo, and the 
system is irreducible [40J. We now prove that such a system cannot be scale- invariant. Suppose 
by way of contradiction that it would be, and pick any fixed p 1. The main theorem in 
[iO] insures that there are two different iable functions a(x) and (3(x) such that the algebraic 
identities: 

a x (x)[f(x) + G{x A )u] + OL y {x)g{x) + a z (x)h(x) = f(a(x),(3(x) 1 x c ) + G(a(x))pu, 
f3 x (x)[f(x) + u] + P y (x)g(x) + /3 z (x)h(x) = g(a(x), /3(x), x c ) 

a(x)a(xc) + /3(x)b(xc) + c(xc) = XAa(xc) + x B b(xc) + c(xc) 

hold for all constant x = {xa,xb,xc) and u, and the vector function x i— )> (a(x), f3(x), z) is 
one-to-one and onto, which implies in particular that 

sup G(a(x)) = G . 
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Dividing by u and taking the limit as u oc in the first identity, we conclude that ol x (x)G(xa) = 
pG{a{x)). Doing the same in the second identity, we conclude that fi x (x) = 0. Finally, taking 
partial derivatives with respect to xa in the third identity: 

a(x c )pG(a(x))/G(x A ) = a x (x)a(x c ) + (3 x (x)b(x c ) = a(x c ) 

is true for all x. Since a(xc) ^ 0, it follows that 

pG{a{x)) = G(xa) 

for all x. We consider two cases: (a) p < 1 and (b) p > 1. Suppose p < 1. Pick any sequence of 
points xW w ith -> G as i -> oc. Then G(a(a;W)) G/p > G, contradicting < G. 

If p > 1, picking a sequence such that G(a(x^)) -> G as z — ^ oc gives the contradiction 
G(x^) pG > G. This shows that the FCD property cannot hold. 
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Supplementary Material 

A characterization of scale invariant responses in enzymatic networks 

5 Circuits that exhibit ASI 

We list here the results of the computational screen as described in the Main Text. After 
showing graphical representations for the 25 identified ASI circuits (21 topologies), we provide 
their equations and parameters. 

For each circuit, four plots are shown: 

(a) a comparison between the plots of XA(t) and #b(£) for the original nonlinear system and 
the respective plots for the linearized approximations, 

(b) the plots showing scale-invariant behavior for step inputs, 

and the comparison between the plots of xc(t) for the original nonlinear system and for 
the quasi-steady state approximation, for 

(c) step input change from 0.3 to 0.36 and 

(d) step input change from 0.5 to 0.6. 
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Circuit 1. 



x A — k uA U 



XB k AB X A 



XC k A CXA 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 



xc 



x c + K AC 



— kp B BXF, 



- k B CXB 



x A 
x A + K BA 

XB 



X B + K FbB 

xc 



xc + K B c 



Parameters: K AB = 0.001191; k AB = 1.466561; K AC = 0.113697; k AC = 1.211993; K BA 
0.001688; k BA = 44.802268; K BC = 0.009891; k B c = 7.239357; K uA = 0.093918; k uA 
11.447219; k AC = 1.211993; K AC = 0.1136927; K Fb = 9.424319; k FB = 22.745736 



Topology 9583[691 1] 



Topology 9583[6911] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 



Output of C node for 1=0.36 
Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 9583[691 1] 




Topology 9583[6911] 



„„,C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 2. 



xa — kuA u 



x B 



x c 



k A B%A 
kACXA 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

XQ 

x c + K AC 



— kp B BXFi 



- k B CXB 



XA _ 

x A + K B a 

XB 



kcAXC 



XA 
X A + K C A 



X B + K FbB 

xc 



xc + K BC 



Parameters: K uA = 0.093918; k uA = 11.447219; K BA = 0.001688; k B A = 44.802268; K C a 
90.209027; k CA = 96.671843; K AB = 0.001191; k AB = 1.466561; K Fb = 9.424319; kp B 
22.745736; K A c = 0.113697; k A c = 1.211993; K BC = 0.009891; k B c = 7.239357 



Topology 8854[6911] 



Topology 8854[6911] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 



'■ Output of C node for 1=0.36 
- Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 8854[6911] 




Topology 8854[6911] 



„„,C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 3. 



7 X A 7 X A 7 X A 

xa k uA u- — — k BA x B kaa%a —r^ — 

x A + K uA x A + K B A x A + K A a 

7 X B %B 7 %B 

xb kabxa- — — ^ kcbXb kbbXb- 



x b + K A b x b + K C b x B + K BB 



7 X C 7 

KBCXB- —p7 KACXA- 



xc + K B c Xq + i^AC 



Parameters: K A a = 7.633962; k A A = 86.238263; K AB = 20.265158; k AB = 5.428752; K A c = 
0.258375; k A c = 62.416585; K BA = 0.003960; k B A = 17.705166; K BB = 31.604578; k BB = 
3.692326; K BC = 44.386408; k BC = 65.027941; K CB = 0.701052; k CB = 26.091557; K uA = 
0.464248; k uA = 1.882348 



Topology 7782[3497] 



Topology 7782[3497] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 7782[3497] 



Topology 7782[3497] 



„„,C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




mmmmi C g (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 4. 



7 X A 7 X A 7 X A 
XA k u AU- — KBAXB k>AA x A 

x A + K uA x A + K B A x A + K A a 

7 X B 1 X B 

xb Kabxa- — — ^ kcbxc — 

x B + K A b x b + K CB 

7 X C 7 X C 

XC KBCXB- KACXA- 



xc + K B c xc + K A c 



Parameters: K A a = 7.633962; k A A = 86.238263; K AB = 20.265158; k AB = 5.428752; K A c = 
0.258375; k A c = 62.416585; K BA = 0.003960; k B A = 17.705166; K BC = 44.386408; k B c = 
65.027941; K CB = 0.701052; k CB = 26.091557; K uA = 0.464248; k uA = 1.882348 



Topology 7863[3497] x -\ 3 Topology 7863[3497] 




Output of C node for 1=0.36 
Output of C node for 1=0.6 



(a) Dynamics of A and B in linearized model (b) Ouput from C nonlinear model 




Time Time 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 5. 



x A — k uA U 



x B 



x c 



k AB X A 
k B CXB 



x A + K uA 

XB 



k BA X B 



x B + K AB 

xc 
xc + K BC 



- kc B xc 



- k AC x A 



x A 
x A + K BA 

X B 



x B + K CB 

xc 
xc + K AC 



Parameters:K A£ > = 63.277600; k AB = 6.638959; K AC = 0.133429; k AC = 55.731406; K BA 
0.011188; k BA = 2.749793; K BC = 0.013374; k BC = 45.175191; K CB = 1.457975; k CB 
2.114949; K uA = 24.589517; k uA = 5.346875 



Topology 7864[1001] 



Topology 7864[1001] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 



8 10 
Time 



(a) Dynamics of A and B in linearized model 




Output of C node for 1=0.36 
Output of C node for 1=0.6 



10 12 14 16 18 20 

Time 



(b) Ouput from C nonlinear model 



Topology 7864[1001] 



C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 




Topology 7864[1001] 




„„,C 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



10 12 14 



(c) Quadratic approx. and output of nonlinear system 



(d) Quadratic approx. and output of nonlinear system 
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Circuit 6. 



7 X A 7 X A 7 X A 7 X A 
X A k uA U ~ —=z k BA X B k AA X A k CA X C 

x A + K uA x A + K BA x A + K AA x A + K CA 

xb k AB x A - — — kcBXc k BB x B - 



x B + K AB x B + K CB x B + K BB 



7 X C 7 

KBCXB- —p7 KACXA- 



xc + K BC XQ + i^AC 



Parameters: K AA = 7.633962; k AA = 86.238263; K AB = 20.265158; k AB = 5.428752; K AC = 
0.258375; k AC = 62.416585; K BA = 0.003960; k BA = 17.705166; K BB = 31.604578; k BB = 
3.692326; K BC = 44.386408; k BC = 65.027941; K CA = 26.714681; k CA = 2.806080; K CB = 
0.701052; k CB = 26.091557; K uA = 0.464248; k uA = 1.882348 



Topology 7053[3497] 



Topology 7053[3497] 




Output from node B/Linearized 

Output from node A/Linearized 

Output from node A Nonlinear 
Output from node B/Nonlinear 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 7053[3497] 



Topology 7053[3497] 



„„,C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




„„,C g (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 7. 



x A — k uA U 



x B 



x c 



k AB X A 
k B CXB 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

xc 
xc + K BC 



- kcBXc 



- k AC x A 



x A 
x A + K BA 
x B 



k AA x A 



x A 
x A + K AA 



kcAXC 



XA 

x A + K CA 



x B + K CB 

xc 
xc + K AC 



Parameters: K AA = 7.633962; k AA = 86.238263; K AB = 20.265158; k AB = 5.428752; K AC 
0.258375; k AC = 62.416585; K BA = 0.003960; k BA = 17.705166; K BC = 44.386408; k BC 
65.027941; K CA = 26.714681; k CA = 2.806080; K CB = 0.701052; k CB = 26.091557; K uA 
0.464248; k uA = 1.882348 



Topology 7134[3497] 



Topology 71 34[3497] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 7134[3497] 



Topology 71 34[3497] 



„„,C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




„„,C g (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 8. 



x A — k uA U 



x B 



x c 



k AB X A 
kACXA 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

XQ 

x c + K AC 



- k FB BX Fi 



- k B cx B 



XA 

x A + K BA 
x B 



x B + K FbB 

xc 



+ kcBXC 



XB 



xc + K BC 



- kccxc 



x B + K CB 

xc 



xc + K C c 



Parameters: K uA = 0.093918; k uA = 11.447219; K BA = 0.001688; k BA = 44.802268; K AB 
0.001191; k AB = 1.466561; K Fb = 9.424319; k FB = 22.745736; K AC = 0.113697; k AC 
1.211993; K BC = 0.009891; k BC = 7.239357; K CB = 30.602013; k CB = 3.811536; K CC 
0.189125; k C c = 17.910182 



Topology 5209[6911] 



Topology 5209[6911] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 



Output of C node for 1=0.36 
Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 5209[6911] 




Topology 5209[6911] 



..... C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 




..... C 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 9. 



xa — kuA u 



x B 



x c 



x A + K uA 



k B AXB 



k A BXA 
kACXA 



XB 



x B + K AB 

XQ 

x c + K AC 



+ kcBXc 



XA 

x A + Kb a 

X B 



kcAXC 



- k B CXB 



XB + K C B 

xc 
xc + K B c 



- k FB BX Fi 



- kccxc 



XA 
X A + K C A 
XB 



X B + K FbB 

xc 



xc + K C c 



Parameters: K uA = 0.093918; k uA = 11.447219; K B a = 0.001688; k B A = 44.802268; K C a 
90.209027; k CA = 96.671843; K AB = 0.001191; k AB = 1.466561; K Fb = 9.424319; k Fe 
22.745736; K Ac = 0.113697; k A c = 1.211993; K BC = 0.009891; k B c = 7.239357; K CB 
30.602013; k CB = 3.811536; K C c = 0.189125; k C c = 17.910182 



Topology 4480[691 1] 



Topology 4480[6911] 




- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 



Output of C node for 1=0.36 
Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 4480[691 1] 




Topology 4480[6911] 



..... C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 




..... C 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 10. 



x A — k uA U 



x B 



x c 



k AB X A 
k B CXB 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

XQ 

xc + K BC 



+ k CB xc 



- k AC x A 



x A 
x A + K BA 

X B 



k AA x A 



XB + K C B 

xc 
x c + K AC 



- k FB BX Fi 



- kccxc 



x A 
x A + K AA 
x B 



x B + K FbB 

xc 



xc + K C c 



Parameters: K AA = 24.989065; k AA = 53.174082; K AB = 0.444375; k AB = 12.053134; K Fb 
1.716920; k Fe = 11.601122; K AC = 0.013988; k AC = 8.521185; K BA = 0.005461; k BA 
7.103952; K BC = 51.850148; k BC = 80.408137; K CB = 5.392001; k CB = 3.086740; K CC 
1.962230; k C c = 17.382010; K uA = 4.387832; k uA = 19.638124 



Topology 5676[305] 



Topology 5676[305] 




- Output from node B/Linearized 

- Output from node A/Linearized 

■ Output from node A Nonlinear 

■ Output from node B/Nonlinear 





Output of C node for 1=0.36 
Output of C node for 1=0.6 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 5676[305] 



Topology 5676[305] 





(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 11. 



x A — k uA U 



x B 



x c 



k AB X A 
k B CXB 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

XQ 

xc + K BC 



+ k CB xc 



- k AC x A 



x A 
x A + K BA 
x B 



XB + K C B 

xc 
x c + K AC 



kp B BXF B 

kccxc 



x B 



x B + K FbB 

xc 



xc + K C c 



Parameters: K AB = 0.444375; k AB = 12.053134; K Fb = 1.716920; kp B = 11.601122; K AC 
0.013988; k AC = 8.521185; K BA = 0.005461; k BA = 7.103952; K BC = 51.850148; k BC 
80.408137; K CB = 5.392001; k CB = 3.086740; K CC = 1.962230; k CC = 17.382010; K uA 
4.387832; k uA = 19.638124 



Topology 5677[305] 



Topology 5677[305] 





- Output from node B/Linearized 

- Output from node A/Linearized 

■ Output from node A Nonlinear 

■ Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 5677[305] 



Topology 5677[305] 




..... C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




..... C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 12. 



XA — kuA u 



X B 



X C 



k AB X A 
kACXA 



XA 
X A + K uA 
X B 



k B A%B 



x B + K AB 
x c + K AC 



+ k CB C 



XA 

x A + K BA 
x B 



+ kcAXC 



- k B cx B 



x B + K CB 

xc 



- k FBB x Fi 



XA 
X A + K C A 

x B 



xc + K BC 



- kccxc 



x B + K FbB 
xc 
xc + K C c 



Parameters: K uA = 0.093918; k uA = 11.447219; K BA = 0.001688; k BA = 44.802268; K C a 
5.026318; k CA = 45.803641; K AB = 0.001191; k AB = 1.466561; K Fb = 9.424319; kp B 
22.745736; K A c = 0.113697; k A c = 1.211993; K BC = 0.009891; k BC = 7.239357; K CB 
30.602013; k CB = 3.811536; K C c = 0.189125; k C c = 17.910182 



Topology 5938[691 1] 



Topology 5938[6911] 



Output from node B/Linearized 
Output from node A/Linearized 
Output from node A Nonlinear 
Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 5938[691 1] 




Topology 5938[6911] 



..... C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 




..... C 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 13. 



x A — k uA U 



x B 



x c 



k AB X A 
k B CXB 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

xc 
xc + K BC 



+ k CB xc 



- k AC x A 



x A 
x A + K BA 

XB 



k AA x A 



XB + K C B 

xc 
xc + K AC 



- k B BXB 



x A 
x A + K AA 

X B 



+ k CA x c 



x A 
x A + K CA 



- k C cx A 



XB + K BB 

xc 
xc + Kcc 



Parameters: K AA = 24.989065; k AA = 53.174082; K AB = 0.444375; k AB = 12.053134; K Fb 
1.716920; kp B = 11.601122; K AC = 0.013988; k AC = 8.521185; K BA = 0.005461; k BA 
7.103952; K BC = 51.850148; k BC = 80.408137; K CB = 5.392001; k CB = 3.086740; K C c 
1.962230; k CC = 17.382010; K uA = 4.387832; k uA = 19.638124; K CA = 15.479253; k CA 
4.903430 



Topology 6405[305] 



Topology 6405[305] 





- Output from node B/Linearized 

- Output from node A/Linearized 

■ Output from node A Nonlinear 

■ Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 6405[305] 



Topology 6405[305] 




,„.,C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 14. 



X A — k uA U 



X B 



X C 



k AB X A 
k B CXB 



X A + K uA 
X B 



k BA X B 



x B + K AB 
xc + K BC 



+ k CB x c 



- k AC x A 



x A 
x A + K BA 
x B 



+ k CA x c 



x B + K CB 

xc 
x c + K AC 



- k BB x B 



x A 
x A + K CA 
x B 



- k C cx A 



x B + K BB 

xc 
xc + K C c 



Parameters: K AB = 0.444375; k AB = 12.053134; K Fb 1.716920; kp B 
0.013988; k AC = 8.521185; K BA = 0.005461; k BA = 7.103952; K BC 
80.408137; K CB = 5.392001; k CB = 3.086740; K C c = 1.962230; k C c 
4.387832; k uA = 19.638124; K CA = 15.479253; k CA = 4.903430 



11.601122; K AC 
51.850148; k BC 
17.382010; K uA 



Topology 6406[305] 



Topology 6406[305] 




- Output from node B/Linearized 

- Output from node A/Linearized 

■ Output from node A Nonlinear 

■ Output from node B/Nonlinear 





Output of C node for 1=0.36 
Output of C node for 1=0.6 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 6406[305] 



Topology 6406[305] 



..... C 3@( t ) b y solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




..... c 6 (t) °y solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 



35 



Circuit 15. 



1 x A x A 
xa k uA u- — — — kbaxb — 

X A + &uA X A + K B A 

7 Xb , X B 

XB = KABXA- KFbBXFb 



x b + K A b x b + K FbB 

, x c 7 x c j x c 

XC k A CXA ~ — —r. k B CXB kccXA — 

x c + K A c x c + K B c x c + K C c 



Parameters: K AB = 0.709169; k AB = 7.445605; K Fb = 1.495375; k Fe = 7.282827; K A c = 
0.002566; k A c = 1.115065; K BA = 0.002522; k B A = 5.753075; K BC = 0.017051; k B c = 
2.777794; K CC = 0.195997; k CC = 1.480130; K uA = 0.225814; k uA = 2.492872 



Topology 3022[5288] x ^ 3 Topology 3022[5288] 




Time Time 



(a) Dynamics of A and B in linearized model (b) Ouput from C nonlinear model 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 16. 



This is the same topology as in the previous case, only a different parameter set was used: 

Parameters: K AB = 0.001191; k AB = 1.466561; K Fg = 9.424319; k Fg = 22.745736; K AC = 
0.113697; k A c = 1.211993; K BA = 0.001688; k B A = 44.802268; K BC = 0.009891; k BC = 
7.239357; K CC = 0.189125; k CC = 17.910182; K uA = 0.093918; k uA = 11.447219 



Topology 3022[691 1 ] x 1 3 Topology 3022[691 1 ] 




Time Time 



(a) Dynamics of A and B in linearized model (b) Ouput from C nonlinear model 



x -| o" 3 Topology 3022[691 1 ] x 1 0" 3 Topology 3022[691 1 ] 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 17. 

This is the same topology as in the previous case, only a different parameter set was used: 

Parameters: K AB = 1.620877; k AB = 2.306216; K Fb = 2.012565; kp B = 2.700847; K AC 
0.010933; k A c = 8.968091; K BA = 0.001812; k B A = 10.039221; K BC = 0.014199; k BC 
17.762333; K CC = 2.686891; k CC = 4.139044; K uA = 0.161715; k uA = 1.933303 



Topology 3022[7239] 



Topology 3022[7239] 



Output from node B/Linearized 
Output from node A/Linearized 
Output from node A Nonlinear 
Output from node B/Nonlinear 





(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 3022[7239] 



Topology 3022[7239] 



..... C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




..... c 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system 
Circuit 18. 



(d) Quadratic approx. and output of nonlinear system 



X A 
±b 
X C 



k u AU 



X A 



k A BXA 
k>AC x A 



x A + K uA 

X B 



k BA x B 



X A 



x B + K AB 

xc 
xc + K AC 



+ k BB x B 



x A + K BA 
x B 



k AA x A 



x A 



k B cx B 



%B + K B b 

xc 
xc + K BC 



kp B BXF 1 



kccxc 



x A + K AA 
x B 



x B + Kf bB 
xc 



xc + Kcc 



Parameters: K AA = 17.569120; k AA = 2.198366; K AB = 9.435176; k AB = 3.134007; K Fb = 
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0.469083; kp B = 1.934194; K AC = 0.062914; k AC = 2.742206; K BA = 0.003245; k BA 
75.352905; K BB = 27.463128; k BB = 10.551155; K BC = 0.041615; k BC = 61.333818; K CC 
0.039332; k C c = 4.756637; K uA = 0.005167; k uA = 8.186533 



Topology 310211212] 



Topology 3102[1212] 



- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 3102[1212] 



Topology 3102[1212] 



..... C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




(c) Quadratic approx. and output of nonlinear system 



(d) Quadratic approx. and output of nonlinear system 
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Circuit 19. 



x A — k uA U 



x B 



x c 



k AB X A 
k B CXB 



XA 

x A + K uA 

XB 



k B AXB 



x B + K AB 

XQ 

xc + K BC 



- k FBB x Fi 



- k AC x A 



x A _ 
x A + K BA 

X B 



k AA x A 



x A 
x A + K AA 



x B + K FbB 

xc 



x c + K AC 



- kccxc 



xc 



xc + K C c 



Parameters: K uA = 4.387832; k uA = 19.638124; K BA = 0.005461; k BA = 7.103952; K AA 
24.989065; k AA = 53.174082; K AB = 0.444375; k AB = 12.053134; K Fb = 1.716920; k Fe 
11.601122; K BC = 51.850148; k BC = 80.408137; K AC = 0.013988; k AC = 8.521185; K CC 
1.962230; k C c = 17.382010 



Topology 3489[305] 



Topology 3489[305] 





- Output from node B/Linearized 

- Output from node A/Linearized 

■ Output from node A Nonlinear 

■ Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 3489[305] 



Topology 3489[305] 




..... C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 



..... C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 20. 



X A — k uA U 



X B 



X C 



k AB X A 
k B CXB 



XA 
X A + K uA 
X B 



k BA X B 



x B + K AB 
xc + K BC 



- k Fj , B x F 



- k AC x A 



x A 
x A + K BA 
x B 



x B + K FbB 

xc 



x c + K AC 



- kccxc 



xc 



xc + K C c 



Parameters: K uA = 4.387832; k uA = 19.638124; K BA = 0.005461; k BA = 7.103952; K AB 
0.444375; k AB = 12.053134; K Fb = 1.716920; k FB = 11.601122; K BC = 51.850148; k BC 
80.408137; K AC = 0.013988; k AC = 8.521185; K CC = 1.962230; k CC = 17.382010 





Topology 3490[305] 



- Output from node B/Linearized 

- Output from node A/Linearized 

■ Output from node A Nonlinear 

■ Output from node B/Nonlinear 



Topology 3490[305] 




- Output of C node for 1=0.36 

- Output of C node for 1=0.6 



0.5 1 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 3490[305] 



Topology 3490[305] 



...„C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 





„„,C 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



1 1.5 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 21. 



X A — k uA U 



X B 



X C 



k AB X A 
kACXA 



XA 
X A + K uA 
X B 



k BA X B 



X A 



X B + K AB 

x c + K AC 



- k Fj , B x F 



- k B cx B 



x A + K BA 
x B 



+ kcAXC 



XA 
X A + K C A 



x B + K FbB 

xc 



xc + K BC 



- kccxc 



xc 



xc + K C c 



Parameters: K uA = 0.093918; k uA = 11.447219; K BA = 0.001688; k BA = 44.802268; K CA 
5.026318; k CA = 45.803641; K AB = 0.001191; k AB = 1.466561; K Fb = 9.424319; k Fe 
22.745736; K AC = 0.113697; k AC = 1.211993; K BC = 0.009891; k BC = 7.239357; K CC 
0.189125; k C c = 17.910182 



Topology 3751 [691 1] 



Topology 3751 [6911] 



Output from node B/Linearized 
Output from node A/Linearized 
Output from node A Nonlinear 
Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 3751 [6911] 




Topology 3751 [6911] 



..... C 36 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 




..... C 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff. eqn. 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 22. 

This is the same topology as in the previous case, only a different parameter set was used: 

Parameters: K AB = 1.620877; k AB = 2.306216; K Fb = 2.012565; k F[3 = 2.700847; K AC = 
0.010933; k AC = 8.968091; K BA = 0.001812; k BA = 10.039221; K BC = 0.014199; k BC = 
17.762333; K CA = 0.002690; k C A = 1.506954; K C c = 2.686891; k C c = 4.139044; K uA = 
0.161715; k uA = 1.933303 



Topology 3751 [7239] 



Topology 3751 [7239] 



Output from node B/Linearized 
Output from node A/Linearized 
Output from node A Nonlinear 
Output from node B/Nonlinear 




(a) Dynamics of A and B in linearized model 



- Output of C node for 1=0.36 

- Output of C node for 1=0.6 




10 12 14 16 18 20 

Time 



(b) Ouput from C nonlinear model 



Topology 3751 [7239] 



Topology 3751 [7239] 



2 4 




C 3 6 (t) by solving quadratic eqn 
C (t) as a solution of nonlinear diff . eqn. 



10 12 
Time 



2 4 



mmmm ,C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




10 12 14 

Time 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 

Circuit 23. 



x A 
x B 

X C 



X A 



k A BXA 
kACXA 



XA + K uA 
X B 



k B AXB 



XA 



X B + K AB 

XQ 
Xc + K A c 



kp BB XF 1 



k B cx B 



XA + K BA 
x B 



kcAXC 



XA 



XA + K C A 



x B + Kf bB 

xc 



xc + K BC 



kccxc 



xc 



xc + K C c 
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Parameters: K uA = 0.093918; k uA = 11.447219; K BA = 0.001688; k BA = 44.802268; K C a 
90.209027; k CA = 96.671843; K AB = 0.001191; k AB = 1.466561; K Fb = 9.424319; k Fs 
22.745736; K AC = 0.113697; k AC = 1.211993; K BC = 0.009891; k BC = 7.239357; K C c 
0.189125; k CC = 17.910182 



Topology 2293[6911] 




Topology 2293[6911] 



- Output from node B/Linearized 

- Output from node A/Linearized 
'■ Output from node A Nonlinear 
'■ Output from node B/Nonlinear 




(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 2293[6911] 




Topology 2293[6911] 



„„,C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 24. 



7 X A 7 X A , 7 X A 

XA KuAU ~ — KBAXB h ^CA^C ~ —f7 

X A + K uA X A + K B A X A + K C A 

/c^B^A- ^7 KFbBXFb 



X B + Kab X B + Kp B B 

/c^C^A- KBCXB- 



X C + ^ AC + K BC 



Parameters: K uA = 0.093918; k uA = 11.447219; K BA = 0.001688; k BA = 44.802268; K C a = 
5.026318; k CA = 45.803641; K AB = 0.001191; k AB = 1.466561; K Fb = 9.424319; kp B = 
22.745736; K A c = 0.113697; k A c = 1.211993; K BC = 0.009891; k BC = 7.239357 



Topology 10312[6911] 



Topology 10312(6911] 



Output from node B/Linearized 
Output from node A/Linearized 
Output from node A Nonlinear 
Output from node B/Nonlinear 





(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 10312(6911] 



...„C 36 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




Topology 10312(6911] 



„„,C 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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Circuit 25. 



This is the same topology as in the previous case, only a different parameter set was used: 

K AB = 1.620877; k AB = 2.306216; K Fb = 2.012565; kp B = 2.700847; K AC = 0.010933; 
k A c = 8.968091; K BA = 0.001812; k B A = 10.039221; K BC = 0.014199; k BC = 17.762333; 
K CA = 0.002690; k CA = 1.506954; K uA = 0.161715; k uA = 1.93330 



Topology 10312[7239] 



Topology 10312[7239] 



Output from node B/Linearized 
Output from node A/Linearized 
Output from node A Nonlinear 
Output from node B/Nonlinear 




Output of C node for 1=0.36 
Output of C node for 1=0.6 




2 4 



10 12 14 

Time 



(a) Dynamics of A and B in linearized model 



(b) Ouput from C nonlinear model 



Topology 10312[7239] 



..... C 3 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff . eqn. 




2 4 



10 12 14 16 18 20 

Time 



Topology 10312(7239] 



..... c 6 (t) by solving quadratic eqn 

C (t) as a solution of nonlinear diff. eqn. 




2 4 6 



10 12 14 16 18 20 

Time 



(c) Quadratic approx. and output of nonlinear system (d) Quadratic approx. and output of nonlinear system 
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6 Ratios XA(i)/xB(i) 

In this section, for each ASI circuit, we show that the ratio XA(t)/xB(t) is approximately 
invariant when inputs are scaled, as discussed in the Main Text. 
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Topology 9583[6911] 



Topology 8854[6911] 



_ratioA 36 /B 3e 
..ratio A B /B K 



ratio A 36 /B 3£ 

6 /B 6 



.36 36 

.....ratio A /B 





Topology 7782[3497] 



Topology 7863[3497] 




Topology 7864[1001] 



Topology 7053[3497] 





Figure 8: XA(t)/xB(t) for circuits 1-6 
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Topology 7134[3497] 



Topology 5209[6911] 




ratio A 36 /B 3£ 

6 /B 6 



.36 36 

.....ratio A /B 




Topology 4480[6911] 





Topology 5676[305] 



ratio A. K /B„ f 

....i ratio A /B 



Topology 5677[305] 



Topology 5938[6911] 





_ratioA ;;c /B ;iei 
....i ratio A 6 /B 6 



Figure 9: XA(t)/xB(t) for circuits 7-12 
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Topology 6405[305] 



Topology 6406[305] 





Topology 3022[5288] 



Topology 3022[6911] 





Topology 3022[7239] 



Topology 3102[1212] 





Figure 10: XA(t)/xB(t) for circuits 13-18 
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Topology 3489[305] 



Topology 3490[305] 




ratio A . /B „ f 

....i ratio A /B 




Topology 3751 [6911] 




Topology 3751 [7239] 




Topology 2293[6911] 





ratio A /B 



Figure 11: XA{t)/xB{t) for circuits 19-24 
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Topology 10312[7239] 

1.3 I 1 1 1 1 1 r- 




\ 1 1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 18 20 

Time 



Figure 12: XA(t)/xB(t) for circuit 25 
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7 Tables 



In this section the following three tables for the 25 identified ASI circuits are shown: 

• Table 1. Relative differences in linearization matrices corresponding to different lineariza- 
tions, Aq.3 = .4.(0.3), .4.0.4 — -4(0.4), . . . , £>o.6 — 13(0-6), rounded to 3 decimal places. The 
corresponding expressions are given by: 



and similarly for £> err . These relative differences are very small. The entries in the table 
are of the following form: A eTT displayed as [an a\2\ a^i CL22} and £> err displayed as \b\ 

• Table 2. Relative error between original (nonlinear) system with an initial state £ = 
(xa, xb) corresponding to u — 0.3, and applied input u — 0.36, and the approximation is 
£ + z(t), where z solves the linear system with an initial condition of zero and a constant 
input of 0.06. Additionally, we provide relative errors between the original (nonlinear) 
system with an initial state corresponding to u = 0.5, and applied input of u — 0.6, 
and the approximation given by £ + z(t), where z solves the linear system with an initial 
condition of zero and a constant input of 0.1. The corresponding expressions are given 



• Table 3. Homogeneity property of the states xa and xb- For a constant input u, it holds 
that a(pu) & pcr(u), where a(u) is a unique steady state (xa,xb)- 



err _ 



w=0.3,0.4,0.5,0.6 
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Circuit 


^err 




1 


[0.069 0.004 


0.005] 


[0.002 0] r 


2 


[0.084 0.006 


0.019 0.015] 


[0.004 0] r 


3 


[0.069 0.004 


0.005] 


[0.002 Of 


4 


[0.114 0.007 


0.011 0.003] 


[0.002 Of 


5 


[0.045 0.003 


0.01 0.033] 


[o op 


6 


[0.075 0.012 


0.021 0.012] 


[0.015 0] T 


7 


[0.057 0.012 


0.021 0.012] 


[0.012 Of 


8 


[0.055 0.012 


0.019 0.009] 


[0.016 Of 


9 


[0.069 0.004 


0.005] 


[0.002 Of 


10 


[0.037 0.022 


0.009 0.0707] 


[0.002 0] r 


1 1 


[0.037 0.022 


0.007 0.009] 


[0.002 Of 


12 


[0.025 0.029 


0.007 0.006] 


[0.012 0] r 


13 


[0.037 0.022 


0.009 0.007] 


[0.002 Of 


14 


[0.036 0.022 


0.007 0.009] 


[0.002 Of 


15 


[0.07 0.004; 0.005] 


[0.002 Of 


16 


[0.07 0.004; 0.005] 


[0.002 Of 


17 


[0.073 0.012; 0.017 0.009] 


[0.015 Of 


18 


[0.051 0.004; 0.005] 


[0.002 0] T 


19 


[0.066 0.013; 0.018 0.009] 


[0.015 0] r 


20 


[0.048 0.013; 0.018 0.009] 


[0.016 Of 


21 


[0.051 0.004; 0.005] 


[0.002 0] r 


22 


[0.233 0; 0.011 0.003] 


[0.002 0] T 


23 


[0.069 0.004; 0.005] 


[0.002 0] r 


24 


[0.051 0.004; 0.005] 


[0.002 0] r 


25 


[0.233 0; 0.011 0.003] 


[0.002 0] r 



Table 1: Relative error in linearization matrices 
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Circuit 


x Amax,u=0.36 


x Bmax,u=0.36 


x Amax,u=0.6 


x B max ,u=0.6 


1 


0.055 


0.011 


0.028 


0.005 


2 


0.008 


0.007 





0.002 


3 


0.055 


0.010 


0.028 


0.005 


4 


0.03 


0.007 


0.012 


0.004 


5 


n Aoi 

0.031 


0.006 


0.003 





6 


0.015 


0.016 


0.011 


0.005 


7 


0.023 


0.021 


0.005 


0.004 


8 


0.023 


0.021 


0.004 


0.004 


9 


0.055 


0.01 


0.028 


0.005 


10 


0.097 


0.020 


0.081 


0.016 


11 


0.010 


0.020 


0.084 


0.016 




n riQQ 
U.Uoo 


U.Uzl 


U.U24 


n m n 
U.U1U 


13 


0.097 


0.020 


0.081 


0.016 


14 


0.010 


0.02 


0.084 


0.016 


15 


0.056 


0.010 


0.028 


0.005 


16 


0.056 


0.010 


0.028 


0.005 


17 


0.027 


0.022 


0.004 


0.004 


18 


0.047 


0.010 


0.028 


0.006 


19 


0.027 


0.023 


0.005 


0.004 


20 


0.023 


0.021 


0.005 


0.004 


21 


0.04 


0.009 


0.034 


0.004 


22 


0.116 


0.027 


0.05 


0.013 


23 


0.055 


0.010 


0.028 


0.005 


24 


0.045 


0.01 


0.027 


0.005 


25 


0.117 


0.03 


0.05 


0.013 



Table 2: Relative error between nonlinear and linearized system 
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Circuit 


cr(u . 3 )/0.3 


cr(n .4)/0.4 


er(iio.5)/0.5 


cr(uo.6)/0.6 


1 


(0.195,0.239) 


(0.193,0.237) 


(0.192,0.236) 


(0.19,0.234) 


2 


(0.199,0.364) 


(0.197,0.359) 


(0.194,0.356) 


(0.192,0.353) 


3 


(0.195,0.239) 


(0.193,0.237) 


(0.192,0.236) 


(0.191,0.234) 


4 


(0.132,0.172) 


(0.131,0.170) 


(0.131,0.169) 


(0.13,0.168) 


5 


(0.591,0.11) 


(0.58,0.109) 


(0.57,0.109) 


(0.561,0.108) 


6 


(0.206,0.526) 


(0.198,0.507) 


(0.192,0.493) 


(0.188,0.481 1 


7 


(0.208,0.529) 


(0.2,0.512) 


(0.194,0.498) 


(0.19,0.486) 


8 


(0.206,0.530) 


(0.199,0.512) 


(0.193,0.499) 


(0.189,0.486) 


9 


(0.195,0.239) 


(0.194,0.237) 


(0.192,0.236) 


(0.190,0.234) 


10 


(0.078,0.083) 


(0.075,0.08) 


(0.073,0.078) 


(0.071,0.076) 


11 


(0.077, 0.083) 


(0.074,0.08) 


(0.072,0.078) 


(0.071,0.076) 


12 


(0.153,0.09) 


(0.145,0.086) 


(0.139,0.082) 


(0.135,0.08) 


13 


(0.078,0.083) 


(0.075,0.08) 


(0.073,0.078) 


(0.071,0.070) 


14 


(0.077,0.083) 


(0.074,0.08) 


(0.072,0.078) 


(0.071,0.076) 


15 


(0.195,0.239) 


(0.193,0.237) 


(0.191,0.235) 


(0.190,0.23-1) 


16 


(0.195,0.239) 


(0.193,0.237) 


(0.191,0.236) 


(0.19,0.234) 


17 


(0.204,0.526) 


(0.197,0.508) 


(0.191,0.494) 


(0.186,0.48) 


18 


(0.196,0.24) 


(0.193,0.238) 


(0.192,0.236) 


(0.19,0.235) 


19 


(0.205,0.528) 


(0.197,0.509) 


(0.192,0.494) 


(0.187,0.481) 


20 


(0.206,0.532) 


(0.199,0.513) 


(0.193,0.5) 


(0.189,0.487) 


21 


(0.196,0.24) 


(0.194,0.237) 


(0.192,0.236) 


(0.191,0.235) 


22 


(0.136,0.177) 


(0.134,0.173) 


(0.133,0.171) 


(0.132,0.17) 


23 


(0.195,0.239) 


(0.193,0.237) 


(0.192,0.236) 


(0.191,0.23-1) 


24 


(0.196,0.240) 


(0.194,0.237) 


(0.192,0.236) 


(0.190,0.235) 


25 


(0.136,0.178) 


(0.134,0.173) 


(0.133,0.171) 


(0.132,0.17) 



Table 3: a(u)/u for constant inputs u = 0.3, 0.4, 0.5, 0.6 
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